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Coherent phase-shift keying (cpsk) and differential phase-shift 
keying (dpsk) are widely used modulation methods in digital com- 
munications. Bandwidth efficiency, good noise immunity, constant 
envelope, and simplicity of implementation make these schemes par- 
ticularly attractive for use over the satellite, terrestrial radio and 
voiceband telephone channels. While system analyses abound in the 
literature, treatment is usually restricted to the additive Gaussian 
channel. Important issues determining ultimate performance, such 
as the joint effect of intersymbol interference and the acquisition of 
carrier phase have not been adequately addressed. The main purpose 
of this paper is to develop analytical tools that can be used to assess 
system performance under practical operating conditions. Pure co- 
herent demodulation schemes such as cpsk are ideals which are 
rarely achieved in practice, and carrier phase must be estimated 
prior to and/ or during data transmission. This requires start-up 
time, as well as added equipment, and the fidelity of the phase 
estimate ultimately determines performance. In contrast, dpsk is 
independent of carrier phase, since decisions are made on phase 
differences. However, this comes at a price, and it is known that ideal 
multiphase dpsk suffers an asymptotic performance penalty of 3 dB 
in signal-to-noise ratio (s/n) over ideal cpsk. We develop a new 
rigorous method for calculating the error rates of both cpsk and 
dpsk, under a variety of operating conditions. In particular, we find 
that the intersymbol interference penalty for quaternary dpsk is about 
1 dB worse in s/n than for cpsk. We demonstrate that the detection 
efficiency of cpsk approaches the ideal, provided that the s/n of the 
phase-recovery circuit is about 10 dB more than that at the receiver 
input. Alternatively, for the same s/n, a 10-baud phase-locked loop 
integration time is required to achieve near-ideal performance. 
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I. INTRODUCTION 

Coherent phase-shift keying (cpsk) and differential phase-shift key- 
ing (dpsk) are two techniques often used in digital communications 
over channels such as satellite, terrestrial radio, and voiceband tele- 
phone. The literature abounds in analyses of their performance under 
a variety of conditions. A sample collection of some of this literature 
may be found in Ref. 1. The chief reasons for the widespread use of 
these techniques are simplicity of implementation, superior perform- 
ance over the additive Gaussian noise channel, minimal bandwidth 
occupancy, and minimal envelope variation. 

The relative performance of cpsk and dpsk systems is well under- 
stood only in the presence of additive Gaussian noise. In this case, the 
detection efficiency of dpsk is known to be about 1 dB (in s/n) below 
that of cpsk for binary modulation and this degradation approaches 3 
dB for multilevel systems. In applications where a 3-dB loss in s/n is 
important, such as in down-link satellite, space communications, and 
terrestrial radio under deep fading conditions, cpsk is the preferred 
method. In cpsk, however, the generation and extraction of a local 
carrier-phase reference at the receiver is required. A coherent phase 
estimate is usually obtained by using phase-locked loop (pll) tech- 
niques, and because of frequency instabilities and phase jitter inherent 
in transmitter and receiver systems, carrier recovery loop bandwidths 
cannot be made arbitrarily small. Consequently, in practice a noisy 
phase estimate is obtained and only partial coherent reception can be 
claimed. The reason for using dpsk is its immunity from slow carrier- 
phase fluctuations; therefore, the phase recovery problem inherent in 
cpsk is avoided. However, the detection efficiency of dpsk may ap- 
proach that of cpsk under noisy phase estimation conditions and 
intersymbol interference (isi). The need to understand this phenome- 
non on a fundamental level is the principal objective of this paper. 

As bandwidth occupancy is always important, the effects of isi 
generated by the use of band-limiting filters must be taken into account 
in any analysis of these systems. Because of the linear nature of the 
demodulation process in cpsk, the effect of isi has been treated in 
great detail. Since dpsk demodulation is inherently nonlinear, the 
analysis of performance is very difficult and no adequate analytical 
methods are currently available. Also, the combined effects of imper- 
fect phase estimation and isi on cpsk must be determined so that the 
relative detection efficiencies of band-limited dpsk and cpsk can be 
fairly assessed. 

In Section II of this paper, we describe a technique for determining 
the degradation in M-ary cpsk operating in the presence of isi, additive 
Gaussian noise, and imperfect carrier phase. In Section III, we consider 
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the performance of M-ary dpsk subject to isi and additive Gaussian 
noise. 



II. COHERENT DETECTION 

2. 1 System description of CPSK 

Figure 1 shows the Af-ary cpsk system that we consider. The signal, 
s(t), before the transmit filter can be represented as 

s(t) = Re{Ax(t)exp[i(27rf c t + /*)]}, i = v^l, (1) 

where the baseband modulation signal is 



x(t)= £ exp(ia k )rect[(t-kT)/Tl 

k 00 



(2) 



and the constants A, f c , and /j are the carrier amplitude, frequency (in 
Hz), and phase, respectively. Also, rect(-) is the rectangular window 
function, T the signaling interval, and the sequence of discrete phases 
[a*] corresponds to the data sequence to be transmitted. Without loss 
of generality, we assume that the M phase values of a* are uniformly 
distributed with equal probability between (— ir, «■]. So, a* takes on 
value in the set a*eA, 
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Fig. 1— M-ary cpsk receiver, M > 2. 
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We also assume that the data phases in different time slots are 
statistically independent. 

In our model, the transmit filter, transmission channel, and receive 
filter are linear and time invariant. Therefore, the complex envelope, 
y(t), at the output of the receive filter may be written as 

y(t)=x(t)®h(t) + n(t) + in(t) 

h(t) = h T {t)G)hc{t)®h R (t), 

where h T {t), hc(t), and h R (t) are, respectively, the impulse response of 
the transmit filter, the channel, and the receive filter. The symbol 
© denotes convolution. Also, n(t) + in (t) is the complex envelope of 
the Gaussian noise passed through the receive filter. For symmetrical 
filters, n(t) andri(J) are independently and identically distributed (iid) 
Gaussian random variables with mean zero, and variance 

o 2 = N i \H R (f)\ 2 df, 



where No is the double-sided spectral density of the original white 
noise and H R (f) is the baseband equivalent transfer function of the 
receive filter. 

2.1.1 Detection in CPSK 

Assuming that the recovered carrier is exp[i(27rf c t + jt )], where p. is 
an estimate of /x in eq. (1), the detector operates on the signal, w(t), 
represented as 

w(t)= J z(t-kT)exp[i(a k + e)] + £+i-n, (3) 

k 00 

where £ and tj are iid gaussian random variables with mean zero, and 
variance a 2 , 

is the phase error, and 

z(t) = h T (t) hc(t) h*(t) rectf ^j . 

To estimate the transmitted phase, a = $eA at t = 0, an ideal cpsk 
detector measures the phase of w(t) at t = to, and a correct decision 
results when 

M M 

9 = phase angle of w(to), 

w(to) = w(t) . (4) 

a -<I> 
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2. 1.2 Error rate for M-ary CPSK 

Here, we briefly review some known results for cpsk and then 
develop new results applicable to our more general model. 

Error-rate calculations for ideal cpsk in added Gaussian noise can 
be found in Refs. 2 to 7. References 8 and 9 provide numerical methods 
for calculating the probability of error in the presence of isi. Reference 
10 takes into account isi and demodulation phase error, but the results 
are restricted to only binary and quaternary systems. We now gener- 
alize these results. 

Using the union bound and the representation of the received signal, 
eq. (3), it follows from eq. (4) that the probability of error, PeflO), 
given that the phase $ is transmitted, is 



max(Pi, P 2 ) < Pe(|$) < P x + P 2> 



where 



Pi = Pr 



sin| - $ + — I < 



= PrJ Im w(to)exp 



-•••- s 



<0 



P 2 = Pr 



sin( 6 - O - — I > 



= Pr^ Im w(to)exp 



-''* + M 



>0 



(5) 



Note that the average symbol probability of error Pe is 

Pe = T? I Peflfc). 

M <t, <A 

But, since the signal constellation is assumed to be circularly symmet- 
ric, Pe(|3>) is independent of 4>. 

For convenience, we shall now assume that <I> = ir/M. Hence, 



Pi = Pr/lmj 



r exp 



.■|A + s + « 



+ £'r*exp 



i ( Bu + ot k + c 



+ T)<0 



(6) 



where 



rifixpdpk) = z(t - kT) 



and Y.' denotes the exclusion of the term k = 0. A similar expression 
can be written for P 2 . 
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Accurate estimations of Pi and P 2 are easy to obtain in the presence 
of only Gaussian noise, but are more difficult when isi is added and 
are even more tedious when the distribution of carrier-phase error, e, 
must be taken into account. 

In the next section we derive an exponentially tight upper bound on 
these quantities for a fixed carrier-phase error and then perform 
asymptotic [large signal-to-noise ratio (s/n)] analyses on these upper 
bounds for a given distribution of carrier-phase error. 

2.1.3 Bounds on the error rate 

We begin by writing eq. (6) as 

> e , (7) 



Pi = < Pr 



1(e) + -q< -rosinf -^ + fa + e j 



where < ) £ denotes expectation with respect to e, and where, 

1(e) = E'r* sin(/3* + <x k + e). (8) 

Before we can proceed with eq. (7), we need specific information on 
the probability density function (pdf) of the demodulating phase error 
e. We shall assume that the phase reference is derived from a pure 
tone by a first-order pll. It is well known 1 that the resulting pdf for 
the phase error, €, is 

exp(X cos e) < . 

p < {€) — ^oiXT' 161 -"' () 

where X is the s/n at the input to the pll multiplied by the reciprocal 
of the pll bandwidth, 

In eq. (10), G is the average power in the carrier, N p is the double- 
sided noise spectral density, and B L is the noise bandwidth of the 
linearized pll. Also, in eq. (9), I (x) represents the modified Bessel 
function of the first kind and of order 0. For a second-order pll, the 
pdf of e is also approximately given by eq. (9). We shall use this density 
to obtain bounds on Pi. 

Since e is a symmetric random variable, eq. (7) yields 

Pl=i<V(€) + V(-€» ei 



where 

1 [rosinUv/M) + fa + e] + 5V r*sin(a* + fa + e) 
Vie) 4 - erfc ^- 
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(11) 



Using upper bounding techniques and Laplace's method," we show in 
Appendix A that 



where 

J\a '■ 



P\ < J\ a + Ji a , 

~ exp(-p 2 [sin 2 [(7r/M) + ft - e ] + D(l - cos €o)]} 
{cos € + (2/D) cos 2[(tt/M) + ft - e ]} 1/2 
exp{-p 2 sin 2 [(7r/M) + ft]} 
{1 + (2/D) cos 2[(ir/M) + ft]} 1/2 



,P 2 = 



rl 



2(o* + oi) 



» 1, 



a 2 I = l'rl,D = - 1 , 
P 



1 • «/ W n 



1 -| COS2 fe + /? °) l 



, D » 1. 



and 



J20 * ( 1 - - ) V2ttZ)p 2 exp[-Z>p 2 (l - cos 8)], S = £ + ft. 

Note that p 2 is the s/n of the system. Also, D can be regarded as the 
ratio of s/n in the phase recovery circuit to that in the psk system or 
the integration time in bauds. 
Similarly, we can show that 

P-z ^ J\a + Jla- 

In summary, the average symbol probability of error, Pe, for M-ary 
cpsk system can be upper bounded by 



Pe< 



exp{-p 2 [sin 2 [(77-/M) + ft - €o ] + D(l - cos e )]} 
{cos eo + (2/D) cos 2[(tt/M) + ft - e]} 1/2 
exp{-p 2 sin 2 [(7r/M) + ft]} 



{l + 2/Z)cos2(77/M + ft)} 1/2 



+ 2 



1- 



(tt/M ) + ft 



sl27rDp- 



x exp< —Dp 2 



I -cos(- + ft 



,p 2 »l, 



1 . J ir n 

€0 = D Sin2 [M + P ° 



l -l co *i^+^ 



Z)» 1. 



This upper bound becomes 
P < 2 exp 



-P 2 sin 2 (- + ft 



(12) 
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when phase estimation is perfect, D -* oo. Equation (12) is the well- 
known Chernoff bound for Af-ary cpsk. 12 
If the observation interval of the pll is large, D » 1, and if M » 1, 



€oS 5 Sm2 fe + /?0 



2 Jtt n 



and 



P < exp -p 2 sin 2 — + fa 



{'- 



M 

2 COS 2 [(7T/M) + fa] 



2 J If n 

1 -D COS2 U + /!0 



+ exp| -p 2 sin 2 ( — + fa 



- exp -p 2 sin 2 — + /?o 



<- 



M 



2 cos 2 [(tt/M) + fa] 
D 



2 J* n 

l -D COS \M + e° 



P 2 = 



r^ 



-, p 2 »l,M>2,Z)»l. 



, ~ -^ -.^^ ■• (13) 

2(a 2 + a 2 ) 

Comparing eqs. (12) and (13), we see that the degradation in s/n 
because of imperfect phase estimate for multiphase cpsk systems is 
asymptotically given by 



G - J1 -! cos2 te + ^ 



where G-* I as D— » oo as it should. 



2 Jit n 



2.2 Example of quaternary (M = 4) CPSK system 

Let us consider a quaternary (M = 4) cpsk system and assume that 
the channel is ideal. 

If 4-pole Butterworth transmit and receive filters are used, the 
resulting average symbol probability of error is plotted in Fig. 2. Note 
that the bound is fairly tight and when the s/n of the phase recovery 
circuit is about 10 dB more than at the receiver input, the detection 
efficiency of cpsk is essentially determined by isi alone. Alternatively, 
we can say that, for the same s/n, a 10-baud pll integration time is 
required to achieve this isi-limited performance. For this filter, the isi 
penalty is about 1 dB. 

If M > 2, it is well known that the penalty in s/n because of Gaussian 
noise alone is asymptotically given by l/[sin 2 (7r/M)]. 
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Fig. 2 — Probability of error for quaternary phase-shift keying (qpsk) with rectangular 
signaling, noisy carrier-phase recovery, and 4-pole Butterworth transmit and receive 
filters. The s/n in decibels is defined as 10 logi [T/2iVo], where No is the double-sided 
noise spectral density and the ideal received signal power has been normalized to unity. 
Parameter D is the ratio of s/n in the phase recovery loop to that in the psk system. 
The double-sided 3-dB bandwidth of the transmit filter is 2/T and that of the receive 
filter is 1.06/T. Sampling time is 1.74T. 

The upper bound in eq. (13) indicates that if the definition of s/n is 
modified to take into account the isi power a 2 , the additional penalty, 
because of imperfect phase estimation, is 



G,= 



sin 2 — + /?o 



M 



x { i -! cos2 fe +/5 °) 



2 Jw „ 

1 -^ C ° S2 M + ^° 



This quantity is plotted in Fig. 3. We observe that the s/n penalty 
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Fig. 3— Signal-to-noise ratio penalty for Af-ary cpsk with imperfect phase estimation. 
Parameter D is the ratio of s/n in the phase recovery loop to that in the psk system. 



because of isi is independent of M. In Fig. 3, also note that G t -» 
l/[sinWM)]as.D->oo. 

III. DIFFERENTIAL DETECTION 

3. 1 System description of DPSK 

The Af-ary dpsk system is shown in Fig. 4. As before, the baseband 
modulated dpsk signal can be represented as 



x(t)= X exp(ia k )rect[(t-kT)/T]. 



(2) 



ft— 



Here, however, the sequence of phases [/?*] = [«* + i - a k ] corresponds 
to the data sequence to be transmitted. Again, we assume that M 
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Fig. 4 — M-axy dpsk system. 

phase values of fS k are equally distributed over the interval [0, 2w) and 
choose 

Pk = (21 - 1) ^, 1 < / < M, modulo 2tt. 

As in cpsk, we represent the set 

IT o7T 7T 

■j-z, — , • • • , (2M - 1) — 
MM M 

by [A]. Also, we shall assume that the phase symbols, fa's, in different 
time slots are statistically iid. 

If the received phasor at time t is indicated by z and the one in the 
succeeding interval is indicated by z d , the detected phase difference 
measured by an ideal differential detector is 

= angle of w, w A z*Zd, 

where * represents the complex conjugate. For the system shown in 
Fig. 4, 



2 - X tgh + ipk)exp(iotk) + n c + in,, 

k 00 



and 



(14) 



z <i = Z (Sk-i + ipk-i)exp(ia k ) + n Cd + in, 

k= — <xi 



where g h and p k are real, 



gk + ipk= | h(t - /x)rect( — - — | dfx. 



(15) 



(16) 
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As before, h(t) is the overall impulse response of the system with 
transfer characteristic 

H(f)=H T (f)H c (f)H R (f). (17) 

In eqs. 14 and 15, n c , n s , n Cd , and n Sd are iid real Gaussian random 
variables with mean zero and variance 



a 2 = No\ \H R {f)\ 2 df, 



where No, is the double-sided spectral density of the added white 
Gaussian noise. In eq. (17), Ht is the transfer function of the transmit 
filter, Hr, of the receive filter, and Hc(f), the transfer function of the 
channel. The assumption that the Gaussian noise at to is independent 
of the noise at to - T can be justified if the receive filter bandwidth is 
small compared with 1/T. Most of our analysis can be extended if 
these two noise samples are correlated. 

3.1.1 Probability of error for M-ary systems 

If the transmitted symbol associated with the time index k = is 
$eA, a correct decision is made when the received phase difference 6 
is such that 

•n _ _ it 
M M 

As before, the following bounds apply: 

max(Pi, P 2 ) < Pe(|$) < Pi + P 2 , 

where 

Pi = Pr 



smi - $ + — I < 



M 



P 2 = Pr 



smi (/ - - — I >0 



M 



(18) 



These statements are identical to the ones that apply to cpsk, but here 
$ represents a "differential phase" and, therefore, the estimation of 
these probabilities becomes extremely involved. Good estimates are 
only available when Gaussian noise is the sole source of impairment. 
We proceed to analyze Pi and observe that a bound on Pi also 
provides a bound on P 2 . Since the calculations are extremely tedious, 
we relegate the details to appendices and strive to develop only the 
main ideas here. Therefore, from eq. (18), we get 
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Pi(|<D) = Pr 


sin^-O+^<0 




= Pix Im w exp 


H-h)] 


..j 


= Pr \ Im z *2rfexp 


i-'(-5); 


<») 


= Pr| 


Re (z)*z d exp 


-i ( O - 


7T 


♦5)] 



<0 



= Pr(Re2f2 2 <0), 



where 



Z\ = z 

22 = 2<*exp 



-'••*-£ + f 



and 2 and z d are given in eqs. (14) and (15). Since 



Re 2*22 = 



Z\ + 22 


2 


2i - 22 


2 


2 



eqs. (19) and (20) yield 
where 



W\ = 



Pi(|3>) = Pr(|u;i|<|ii;2|), 
2i + 2 2 2 + 2 rf exp - i[$ - (tt/M ) + (w/2)] 



= X - i (g* + ipk) + (gk-\ + ipk-\) 



X exp 



... 77 77 



exp(ia*) + 1+ + *V 



a>2 = 



2i - 2 2 2 - 2,/exp - i[0 - (it/M ) + (tt/2)] 






x exp 



., , IT IT 

-•'*-S + 2 



exp(i'a A ) + |- + it}-, 



(19) 



(20) 



(21) 



and where £+, 7/+, £_, and tt_ are Gaussian noise terms, given by 
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£+ = 



n c + Re{n C d + in S d) exp - i[0> - (it/M) + (tt/2)] 



y+ = 



n s + Im(n C rf + in sd ) exp - i[$ - (tt/M) + (tr/2)] 



£- = 



ric - Re(n cd + in S d) exp - i[4> - (tt/M) + (tt/2)] 



and 



V- 



n s - Im(n crf + triad) exp - i[$ - (tt/M) + (it/2)] 



It can be verified that the above Gaussian random variables are iid 

with mean zero and variance o 2 /2. 

3. 1.2 Exact computation of probability of error 

For a given symbol sequence, the conditional probability of error is 
seen from eq. (21) to be given by the probability that a particular 
Gaussian quadratic form exceeds another. This is a well-known prob- 
lem and the answer can conveniently be expressed in terms of the 
tabulated Marcum Q function. 13 Thus, after some algebra, eq. (21) can 
be shown as 14 



Pi (|<£, symbol sequence) 



or 



o? + (A 



l-Q 



a- 



+ ^Uq 



vaf + a| va? + <r|/ . 
a- a+ 



where 



oi + ffi VVa? + <ri ' Va? + <&/ 



a + x' 
Q(a, b) = I expl )I (ax)xdx, 



(22) 



and /„(•) is the modified Bessel function of the first kind and of 
order n, 

a+ = \(wi)^ + \, 
a- = |<ufe)e_,v_|i 
of - <£> - <*> - «*/2 



and 



oi = <£> = (V-) = o*/2. 
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The major difficulty at this point is clearly carrying out the averages 
in eq. (22) over all possible symbol sequences. In general, a- and a+ 
contain an infinite number of isi terms and the averaging process is 
difficult. Clearly, for a small number of isi terms, it can be carried out 
by enumeration. But, in general, the number of terms in computing 
the average explodes exponentially and enumeration becomes intrac- 
table. For example, for 10 isi terms and a quaternary dpsk system, the 
number of terms is about a million! So, we obviously need more 
efficient methods of estimating these averages. 15 

In this paper, we assume that the number of dominant isi terms 
contained in a+ and a_ is not large and that they become insignificant 
when isi samples are far away from the desired sample. Assuming that 
the same number N indicates the number of dominant preceding and 
succeeding isi samples (total significant isi terms is 2N), our approach, 
then, is to obtain upper and lower bounds on Pi as a function of N and 
demonstrate that these bounds coincide with N — > <». 

For any N, the evaluation of these bounds requires M 2N computa- 
tions. This can be carried out with modest effort on a high-speed 
digital computer. The error becomes smaller when N is increased. 

We show in Appendix B that the error probability can be bounded 
as 



X 1 (N)<P 1 (|0>)< X 2(iV), 



where 



* W = 



+ 



1 + (1 - A) 2 
(1 - A) 2 



v/2 



a+ 



\/2a-(l - A) 



_Wl + (1 - A) 2 ) ' aVl + (1 - A) 2 . 
\/2a-(l - A) V2a 



1 + (1-A) 2 \ _ ff Vl + (l-A) 2 'o>/l + (l-A) 2 _ 



1- (Q 



— 4 expl — 



•ha- >/2p 



AV 



— 4 expl — 



Ay 



2 I Hi 

k<-N 
k>N 



2 2 Gi 

k<-N 



(23) 



and 
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X2(N) = 



1 + (1 + A) 2 
(1 + A) 5 



V2 



a+ 



>/2a-(l + A) 



.ay/1 + (1 + A) 2 ' aVl + (1 + A) 2 . 
V2a-(1 + A) V2a_ 



+ 1 + (1 + A) 2 \f . Wl + (1 + A) 2 ' aVl + (1 + A) 2 /. 



1- (Q 



V2a- V2p 



i a 



+ 4 expl — 



AV 



2 2 Gl 



+ 4 expl — 



AV 



2 2 Hi 

k<-N 



(24) 



In eqs. (23) and (24), A andp are arbitrary, < A < 1, and 



-il* 



G* = -](^*+ ijp*) + (#a-i + ip*-i)exp 



-IH 



r 1 (gk + ipk) - (gk-i + ipk-i)exp 



-M + 2J_j 



For any TV, we can choose A and p by trial and error so that the 
difference between the upper and lower bounds is a minimum. Since 
this optimization is not critical to our method, we choose 



P = 



2a 2 */<a 2 -> . o 2 . 
A V a 2 or 



-1 1/2 



where 



or = max( 2 Gl I Hi). 

\k<-N k<-N / 



k>N k>N 

For A «: 1, the difference, Z, between the upper and the lower bounds 
can be shown as 



Z=2A 



1 + 



a: 



a: 



a 2 o 2 

-la 2 -)\ o% ({a 2 -) 



Pi(\®,N) 



+ 8exp {-T^fK^ +ln 5 A )]} 



+ 8 exp\ - 



k<-N 

A 



2 Hi 

k<-N 
k>N 



oll-^ + ln^-A 



(25) 
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where 



Pi(|*,W) 
1 
~2 



1 "<«T-T 



and a+ and a_ contain only the first 2iV significant terms. When N —* 
oo, Z in eq. (25) can be seen to approach zero. 

Since a+ and a- in eqs. (23) and (24) contain a finite number of isi 
terms, we can use the direct method to evaluate the averages and then 
compute the bounds. We choose the initial N so that Or < 1, A = 
sor. We then increase N so that the desired accuracy of computation 
is achieved. 

3.1.3 Upper bound on the probability of error 

Since the exact evaluation of Pi(|4>) is difficult — though we have 
developed in the last section numerical techniques which can be used 
to compute Pi with any desired accuracy — we attempt to derive an 
upper bound on Pi. 

Although our bounding approach seems reasonable, the final bound 
that we obtain turns out to be loose. Our purpose in including this 
section is to alert readers about this approach and to emphasize the 
importance of the tedious, but necessary, computations outlined in 
Section 3.1.2. 

To facilitate our bounding techniques, we need the following rela- 
tions. For any two random variables x and y and any two real numbers 
a and A, we can show (see Fig. 5) that 

Pr(* > a + A) - Pr(;y < -A) 

< Pr(x + y > a) 

< Pr(x > a - A) + Pr(;y > A). (26) 
Equations (21) and (26) yield 

Pi(|O) = Pr(|iife|-|a;i|>0) 

< Pr(| W2 1 > A) + Pr(-| uh \ > -A) 

= Prfliifel > A) + Pr(|i£>i| < A), (27) 

where A is arbitrary. We choose A > so that the upper bound in eq. 
(27) is a minimum. The method of choosing A will be discussed later. 
Now, for any complex random variable z = x + iy, we can show (see 
Fig. 6) that 
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Fig. 5 — (a) Upper bound on Pr(* + y > a), x and y, any two random variables and A 
is arbitrary, (b) Lower bound on Pr(.x + y > a), x and y, any two random variables and 
A is arbitrary. 

Pr(|z|>a) 

<Pr(|Re2|>ai) + Pr(|Imz|> Va 2 - a?). 
Hence, 
Prfliifel >A) 

< Pr(| Re w 2 1 > Ai) + Pr(| Im w 2 \ > A 2 ), A\ + Al = A 2 , (28) 
where 

Re w 2 = £- + J CkCos(a k + A*) 



Im w 2 = v- + S Cksin(a k + \k) 



k 



Ck= S 1 Cft + *P*) _ (£*-i + »PA-i)exp 
and 
CAexpd'A*) = - | (^* + ipk) - (gk-\ + ip*-i)exp 



-i $ 



7T 7r\~| I 

-M + 2J.j 

77 7r\ I 

M + 2J.} 



-a # - s+ . 



Since the real and imaginary parts of w 2 are the sum of a Gaussian 
random variable and a set of interference terms, various methods given 
in Refs. 16 to 18 and 19 to 27 can be used to bound Pr(| Re w 2 \ > A\) 
and Pr(|Im w 2 \ > A 2 ), A 2 + A 2 = A 2 . Even though the other bounds 
are sometimes claimed to be tighter, we shall use the simpler Chernoff 
bounds. 
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Fig. 6— Upper bound on Pr(| z\ > a). Parameter A satisfies < A < a. 



Appendix C shows that 
Pr(|u> 2 | >A) < 2 exp 

where 



(A - Cmi) 5 
o^ + o 2 - 



+ 2 exp 



(A 2 — CW' 



<? + <£ 



( TT IT 

M=A C ° S (4 + TM 

. / 7T 77 \ 

Ai = A sin - + — — 



v 4 2MJ 
Cmi = max{[C cos(a + Ao) + CiCOs(ai + Xi)]} 



and 



ai - cto = — 
M 



Cm2 = max{[C sin(ao + Xo) + Cisin(ai + Ai)]} 



ai — ao = — 
M 



a 2 - = y:'c\ 



A=-oo 



Also, for any complex random variable z = x + iy, we can show (see 
Fig. 7) that 
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Hence, 



Since 



and 



— |lm/|<3 




Im; 







|Re z|< a 



Fig. 7 — Upper bound on Pr(| z\ < a). 

Pr(\z\<a) <Pr(|Re2|<a, |Im z\ < a) 
Pr(|z|<a) <Pr(|Rez|<a) 
Pr(|z|<a) < Pr(| Im z \ < a). 

Pr(| w x | < A) < Pr(| Re u>i | < A). 

Re w\ < | Re W\ \, 
Pr(| Re wi | < A) < Pr(Re u>i < A), 

Pr(| h>i I < A) < Pr(Re u>i < A). 



We write 

Re Wi = £+ + X D k cos(ctk + 5a) 



Im wi = j]+ + £ -DAsin(a A + S k ), 

k=-oo 



D k = 



- \ (g k + ipk) + igk-i + ipk-i)exp 



-ilQ 



D k exp(i8 k ) = - | (g k + ip*) + (gk-i + *p*-i)exp ~M ° ~ ^ + f ) |> 
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and 



£ D k cos(ctk + 8 k ) 



= D cos(a + 8q) + Acos(ai + Si) + £" D k cos(a k + S k ). 
Using Chernoff bounding techniques, it can be shown that 

(Dm - A) 2 ' 



Pr(| wi | < A) < exp 



<? + <£ 



where 



Dm = min[DoCOs(ao + So) + Acos(ai + Si)] 

«1 — CKO = rr-L 

M 
The upper bound on Pi(|/3o = ir/M) can, therefore, be written as 



Pi 



A> = ^)<2exp 



(A,- 


- Cmi) 






^H-o 2 - 




+ 2 exp 


(A 2 

a 2 


- CjVf2) 


+ exp 



<? + <£ 



Ai - C m i > 0, A 2 - C M 2 > 0, Dm - A > 0, A\ + A\ = A 2 . (29) 

The bound is minimum when A, A\, and Ai are chosen so that the 
derivative of eq. (29) is zero. This can be found by using well-known 
numerical methods. 

3.2 Example of quaternary (M = 4) DPSK system 

Let us consider a quaternary (M = 4) dpsk system and assume that 
the channel is ideal. 

If 4-pole Butterworth transmit and receive filters are used, the 
bound given by eq. (29) is plotted in Fig. 8. The bound with zero isi is 
plotted in Fig. 9. The exact probability of error with isi is plotted in 
Fig. 10. With or without isi, the bound is unfortunately not very tight. 
Actually, one can show that the penalty as predicted by the bound 
with zero isi is about 4.6 dB worse than the actual penalty for a binary 
system, and 8.3 dB worse for a quaternary system. This is inherent in 
our techniques and not the result of using Chernoff bounding methods. 
In our opinion, obtaining tighter bounds is still an open problem. 
Comparing Figs. 2 and 10, we note that isi penalty for quaternary 
dpsk is about 1 dB worse than for cpsk. We needed 9 isi terms to 
compute Pe with 5 percent accuracy. 
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I I I I I I I I I 



UPPER BOUND FOR DQPSK 



-WITH ISI 




22 



24 



26 28 30 32 

SIGNAL-TO-NOISE RATIO IN DECIBELS 



34 



36 



Fig. 8 — Upper bound on the probability of error for differential qpsk (dqpsk) with 
the same transmit and receive filters as in Fig. 2. Other assumptions are as in Fig. 2. 



IV. SUMMARY AND CONCLUSIONS 

For multiphase M-ary cpsk, we develop an analytical procedure for 
determining detection efficiency when the system is subject to additive 
Gaussian noise, isi, and imperfect carrier-phase estimation. For a large 
s/n, we provide a simple formula for calculating the combined penalty 
caused by isi and noisy phase recovery. For multiphase dpsk, where 
the detection is inherently nonlinear, a rigorous method is developed 
for calculating the error rate in the presence of isi and additive 
Gaussian noise. Using these analytical techniques, it is possible to 
compare the performance of cpsk and dpsk and examine various 
parameter trade-offs. Numerical examples are provided to illustrate 
our methods. 
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10 



"i i i r 



10- 



~l — ' — \~t~2 

DQPSK WITH ZERO 
ISI 




UPPER BOUND- 



- EXACT 



I 1 L 



14 16 18 20 22 24 26 2B 

SIGNAL TO-NOISE RATIO IN DECIBELS 

Fig. 9— Upper bound on the probability of error for DQPSK with zero ISI. 



APPENDIX A 

Chernoff Bound on the Probability of Error 
From Section 2.1.3, 



P, =i | |V(f) ! V(-c)\p.U)rU 



where 



[V(e)+ V(-€)]p,(e)de 
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13 15 17 19 21 23 

SIGNAL-TO-NOISE RATIO IN DECIBELS 



25 



Fig. 10 — Probability of error for dqpsk with the same transmit and receive filters as 
in Fig. 2. Other assumptions are as in Fig. 2. Note that 9 isi terms were needed to 
compute Pe with 5 percent accuracy. 



-r 

-r 

Js 



/ - I [V(e) + V(-e)-\p t (e)de 
[V(<0 + V(-e)]p e (e)d€, 



(30) 



and where V(e) is given in eq. (11). Note that sin{ir/M + /? - e) > for 
< e < tt/M + j8 ; also, sin(?r/M + /? + e) > for < e < <n - (tt/M 
+ Po). Hence, sin(7r/M + /? + e) > for < e < 8, where 



8-min[£+/fe, *-(£ + *)]. 
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Since 



< erfcU) < 2, 

0<[V(e) + V(-e)]<2, 
and, therefore, 

J 2 <2 ) p £ (e)de = Pr(|e|>5). 
Now, from eq. (11), 
V(e) = Pr 



-tj - 7(e) > r sin( — + /Jo + e 



(31) 



where tj is a zero mean Gaussian random variable with variance a 2 and 
7(e) is given in eq. (8). 
Using the Chernoff bound 



Pr[* > a] < exp(-/ia) < exp(jix) ) , /i > 0, 
eq. (31) yields 

V(e) < exp 



/it 
— Ar sin -77 + /fo + € 
\M 



exp{-A[Tj + 7(e)]}. (32) 



For a given e, 7 and tj are independent, and since the data phases 
<Xk in different time slots are iid, 

\2 2 

exp{-A[Tj + 7(e)]} = exp — - \[ < exp[-Ar*sin(/?* + a k + e)]> „. (33) 

We shall now assume that M is an even number so that if <I>eA, 
(77 + 4>)eA. Hence, 

(exp[-Ar*sin(/?* + a k + e)]) ttlt , < a* < 2tt 

= ^ (exp[-Ar*sin08* + a k + e)] 

+ exp[-Ar*sin(/?* + 77 + a* + e)])„ A , < a k < 77 
= (cosh Ar*sin(/?* + a* + e))a A > < a* < 77. 



Since 



cosh x < exp(x 2 /2), 

(exp[— Xrksmiflk + a k + e)])„ A , < a* < 277 



< (exp 

< exp 



A2..2 

A rk . 



2 
A rk 



sin (/?* + ak + e) 



>„., < a* < 77 



(34) 
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From eqs. (32), (33), and (34), 
V(e) < exp 
where 



-Ar sin — + fa + e 



exp 



jtf + rf) 



A>0, 



Similarly, 

V(-e) < exp 



a 2 / = S'ri 



lit 
— Ar sin( — + fa — e 
\M 



exp 



j(o 2 + o$) 



Hence, for < e < tt/M + fa, sin(w/M + fa - e) > 0, and 
xn f rfein 2 [(7T/M)+fa-e] 

V{ - e) " exp WT3F) 

Also, for < e < it - (tt/M + fa), sin(vr/M + fa + e) > 0, and 

rfam 2 [(7r/M) + fa + e] 



V(e) < exp 



Hence, 

V(e) + V(-e) < exp 



2(o 2 + of) 

rfcm 2 [(n/M) + fa + e] 



+ exp 



2(<? + a?) 
rgsin 2 [(7r/M) + /?o - e] 
2(a 2 + of) 



= exp 



-p 2 sin 2 ( — + fa + e 



+ exp 



P 2 = 



t- 2 

^0 



2(a 2 + a?) ' 
< € < 8 = min 



M 



_ ps i n *(_ + /3 - 6 
a? = £' ii 






(35) 



The parameter p 2 is the s/n of the system. 

Note that fa is the phase angle of the complex overall impulse 
response evaluated at t = to. In a well-designed system, it is usually 
small, and eq. (35) shows that the optimum value of fa is zero. Also, 
since fa is usually small and we are interested in M > 2, 6* is usually 
■n/M + fa. 

Thus, from eqs. (9), (30), and (35), we conclude that 
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Jl < J la = 



1 



2irh(Dp*) 



-p L sin"| — + /?o - e 



+ exp 



exp 



-p 2 sin 2 ( — + yS + e 



exp[Z)p 2 cos e]de, (36) 



where the quantity D can be regarded as the ratio of s/n in the phase 
recovery circuit to that in the psk system or the integration time in 
bauds. 

We have not been able to evaluate eq. (31) in closed form but, if 
desired, numerical techniques can be used. To obtain physical insight, 
we shall assume that p 2 » 1 and use Laplace's method to evaluate eq. 
(36); the technique is an application of the following theorem: If h(t) 
is a real function of a real variable t, has a unique maximum at t = 
a, ai < a < ot2, and if x is a large positive variable, it can be shown 
that 11 



fix) = g(t)exp[xh(t)]dt as g(a)exp[xh(a)] 



—n 



2xh"(a) 



1/2 



From 


eq. 


36), 






•J\a 




1 




2ttIo(Dp 2 ) 




J b 


■f 

Jo 


expj p 2 




Jc 


-f 

Jo 


expj p 2 



{J,, + Je), 



D cos e — sin 2 ( e — — — /?o 



D cos e - sin e + — + B 

M 



de, 



de. 



(37) 



(38) 



(39) 



The saddle point € at which the exponent in eq. (38) reaches its 
maximum in (0, 5) is given by the solution of 

1 



sin € 



sin 2[(tt/M) + /?o - eo] D 



(40) 



The transcendental eq. (39) can only be solved numerically. However, 
we can obtain a series solution for eo by using Lagrange's reversion 
formula. 27 If a function f(z) is regular in a neighborhood of z Q and if 
f(zo) = Wo, f'(zo) ¥> 0, then it can be shown 28 that 

f(z) = w 
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has a unique solution 



..^jJSz-s^™,.} , 



where 



n =i ni 



*(z)« 



(41) 



*-'0 



Z — Zq 

f(z) - Wo 



Choosing Zo = 0, eqs. (40) and (41) yield 
1 • J* n 

€o= £ sin2 U +/?o 

From Laplace's formula and eq. (38) 



D » 1, < e < 5. 



J b « exp{ -p 



sin — + fa — €o 



— Z) cos € \ 

2p\D cos eo + 2 cos 2(tt/M + fa - *>)] J " 

Similarly, it can be shown that the exponent in eq. (40) reaches its 
maximum in (0, 5) at e = and 



J c « expl -p 



sin 2 — + fa - Z) 



M 



1 1/2 



2p z [Z) + 2 cos(w/M + #>)] 



For p 2 » 1, D > 1, 



/o(iV) a 



^expCDp 2 ) 



V2ttZ)p 2 ' 

From eqs. (37) to (39) and (42) to (44) 

^ 1 / exp{-p 2 [sin 2 [(7r/M) + fa - e ] - D(\ - cos e )]} 
loS 2\ [cos eo + (2/D)cos 2[(tt/M) + fa- e ]] 1/2 

exp[-p 2 sin 2 |>/M ) + #>)] 



(43) 



(44) 



+ 



[1 + (2//>)coa 2[(*r/M) + /8b]] 



1/2 J" 



Also, 
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J 2 < J 2a = 2 p t (e)dt 
h 



I pA 

J6 



7T/o(Z>p 2 )_ 

•rr — 8 exp(Z)p 2 cos 5) 
~~ ff / (Z>p 2 ) 



exp(Z)p 2 cos erfe) 



= ( 1 - - j V27rDp 2 exp[-Dp 2 (l - cos 5)] 



p 2 » 1, Z> > 1. 



APPENDIX B 

Upper and Lower Bounds on the Probability of Error in DPSK 
Let us write 

z = zn + n c + in s + Zr 
Zd = ZdN + n Cd + in Bd + ZdR, 



where 



and 



zn = £ (gk + ipk)exp{icc k ) 

kz-N 

k^N 

zr = J (gk + ipk)exp(ia k ) 

k<-N 
k>N 

ZdN = £ (£*-i + ipk-i)exp(ictk) 

k>-N 
k^N 



ZdR = £ igk-i + ipk-i)exp{ia k ). 

k<-N 
k>N 



Note that zn and ZdN contain a finite number of isi terms, whereas 2/? 
and ZdR contain an infinite number. Without loss of generality, we shall 
also assume that gk and pk are monotonic decreasing functions of | k |. 
N is an arbitrary positive number. 
Now 

Pi (| $) = Pr(| win + Wut | < | w 2 N + w 2 r |), (45) 

where 
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zn + 2rfivexp 






Win = 



WlR = 



W 2 N = 



and 




+ £+ + hj+, 



+ £- + i-q- 



W2R = 

Note that if in and W2N contain a finite number of isi terms, and £+ + 
gn+ and £- + irj- are independent zero means complex Gaussian 
processes. 
For any two complex numbers si and s 2 , 

|si| - |s 2 |< |si + s 2 |< I si I + \s 2 \. (46) 

Hence, eqs. (45) and (46) yield 
Pr(\w 1N \<\w 2N \ - | M>i*| - \w 2 r\) <Pi(|$) 

= Pr(| Win + w iR \ < \ w 2 n + w 2 r |) 

< Pr(|u;iiv| < \w 2N \ + \w\r\ + \w 2R \) 

or 
w j ! M <1 _ | U » UI |-t-|u to L 



W 2N \ 



W2N\ 



< r' l .(lf^d< i , 1^1 + 1^1 



w 2N \ 



Using the bound in eq. (26), eq. (46) yields 

Pr | jifH< 1 - aV Pr f l "" fl| + I "' 2 '' l > A I < /'."•, 



\W2N\ 



W 2 n\ j \ \W 2 N\ 

£Pr (i^< 1+4 ) +Pr fi2id±l«siL>A ,47) 

where we choose < A < 1. 
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For any two random variables x and v and any a > 0, 10 

Pr(xy>a) <PrMx| >-) + Pr(\y\ >p), (48) 

where p > is arbitrary. From eqs. (47) and (48), 
p r (ii^i< i - a) - Pr(\w 2N \ <p) 

- Pt(\wir\ + \w 2 r\ > Ap) <Pi(|0>) 
< p r (|i^< i + A) + Prd^l <p) 

\\W2n\ J 

+ Pr(|u; li? | + |u; 2 *|>A/>). (49) 

Since vom and wzn are independent complex Gaussian processes, 
from Ref. 13, 



p r |l^< 1± A) = 



U)2N\ 



oi 



<Q 



oi + old ± A) 2 
a+ 



a-(l ± A) 



. Vo? + a 2 2 (l ± A) 2 ' Va? + aid ± A) 2 . 
ai(l ± A) 2 



>h 



+ 



x<Q 



a? + aid ± A) 2 
a-(l ± A) a. 



„2_„2_° 

>,oi-o 2 - — , 



(50) 



.Va 2 + ai(l±A) 2 ' Va? + ai(l±A) 2 . 

where a+ and a- are the appropriate truncated values of a+ and a 
defined in Section 3.1.2. Also, from Ref. 13, 



13, 

Pr(\w 2N \<p) = l-<Q (—,?-)>. 

\ a 2 a 2 / 

Using Chernoff bounding techniques, we can also show that 

/ A 2 p 2 

Prd^ > Ap) < 4 expf- 2 ^ 

k<-N 
k>N 

Pr(|u*«| > Ap) < 4 ex p(~ 2 f^fff) ' 

k<-N 



(51) 



(52) 



(53) 



k>N 



where 



- j (gk + ip*) + (^a-i + jp*-i)exp -i I $ - ~ + ^ J I 
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- \ (gk + ipk) - (gk-i + ip*-i)exp 



and 
H k = 

Equations (49) to (53) yield 

X i(N)<Pi(|0)<X2(N), 
where 



-»!• 



TT 7r\ 1 



*' (iV) = 1 + (1 - A)' 



{- 



(1 - A) 2 
1 + (1 - A) 



A 



a+ 



x/2a_(l - A) 



*> 



.ay/1 + (1 - A) 2 ' aVl + (1 

%/2a-(l - A) >/2a+ 

LaVl + (1 - A) 2 ' aVl + U-A) 2 ^ 

, V2a- >/2p N 
1- <QI ,— -!> 



<Q 



a a 

k 2„2 



— 4 exp / — 



Ay 



2 I Gi 

k<-N 



- 4 exp / — 




and 

X2(N) = 



L + (1 + A) 2 



«>- 



A 



a+ 



y/2a-{l + A) 



.ay/1 + (1 + A) 2 ' aVl + (1 



ttfM 



(1 + A) 2 I " V2a-(1 + A) \/2a + ~ | 

+ 1 + (1 + A) 2 <Q .aVl + d + A^'oVl + d + A) 2 ^ 

,'>/2a- V2p 
+ U-<Q| .— 1> 



a o 

l2„2 



+ 4 exp / — 



Ay 



2 £ G* 

A<-iV 
k>N 



+ 4 exp [ — 




where A and p are arbitrary. 
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APPENDIX C 

Upper Bound on the Probability of Error 

Now, 

Pr(| Re w 2 1 > Ai) = Pr(Re w 2 > Ai) + Pr(-Re w 2 > Ax) 
and 
Pr(Re w 2 > Ai) < exp(-fiAi) <exp(/i Re w 2 )> 

= exp(-/iAi)<exp(/n ite w 2 )>. 

Since the Gaussian random variable £_ is assumed to be independent 
of the interference 



<exp(/x Re w 2 )) = exp(/i 2 a 2 /4) exp 



ju £ C*cos(a A + X*) 



Since /?* = a*+i — a*, and we assume that /fo's are iid, we can assume 
that oik's are independent, and 



a* £ A, k even 



a* 



„ 2tt 4n it 

MM M 



A A s , k odd; 



that is, the signal constellation for odd (even) k can be obtained by a 
simple rotation of the constellation for even (odd) k. 

Let us now assume that the transmitted symbol 3> is tt/M so that 



,>,(l<i>) = />,( 1ft. JjL 



-»2*(|*-«.-b)- 



Noting that 



£ C*cos(a* + A*) = C cos(a + X ) + CiCOs(ai + Xi) 



k — 



+ X" C*cos(o* + X*), 
Pr( Re w 2 > Ax | otx - a = — I < expj -fiAi + ^- 

+ ju[C cos(a + Xo) + CiCos(ai + Xi)] \ 

X (exp[/x X" C*cos(a* + X A )]>. 

We use the notation " to indicate that k = 0, and k = 1 terms are not 
included. 
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Now, 

<exp[*i £" C k cos(a k + A*)]> = II" < exp[/iC*cos(a* + A*)]> 

= n" <exp[/iC*cos(a* + \ A )]> 

even 

X II" (exp[[iCkCos{a k + A*)]). 

odd 

Most often, M = 2 L , L an integer, and since this assumption simplifies 
our bound, we shall assume that M is an integer power of two. (A 
slightly more complex bound can be derived if M ¥> 2 L .) Now, 

H" (exp[/iGfcCOs(a* + A*) ])k !-*«»*==*) 

even 

= H" (COSh[jUC*COS(a* + A*)]>(a*|0<aA<i7) 



= U" - <cosh[/iC*cos(a* + A*)] 

pvpn — 



+ cosh 



IT 



juCaCOs) - + ak + A* 



> 



(*H"-s) 



= n" < c ° sh 

even 

X cosh 



•R I TT 

fiC kcos — cos I ak + Ak + — 



_ . IT . I - IT 

/iC*sin— sin I ak + A* + — 



- («*l<x«sj) 



Since 



and 



<cosh 



cosh x < expfljej), 



cosh x < exp(x 2 /2), 



[iCkCOS — cosl ak + Aa + — 



X cosh 



liCkCos — sinl a* + A* + - 



> 



< exp(/iC*) 



and 

<cosh 



fiC kcos - cosl a* + A* + - 



X cosh 



_, 77 . / . TT 

juCacos- sml a* + A* + — 



> 



< exp 






(u»|0*a»*w/2) 
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Identical bounds can be derived when k is odd. Therefore, we have 
<apD£» C*eos(a* + X*)]> s exp ( /i £"C* + ^- Y>" Cl\, 

\ Aefi, 4 A e nc y 

where ft i is a subset of [••■ , —2, — 1, 2, 3, •••]. For simplicity, we 
choose Qi\ to be the null set. 



Choosing optimum /i, 



Prl Re w 2 > Ai\cti - cto = — 



< exp 



{Ai - [C cos(a + A ) + CiCOs(ai + Ai)]} : 



{a 2 + r CI} 
Ai - [Cocos(a + Ao) + CiCos(ai + Ai)] > 0. 



Similarly, we can show that 
Pr( -Re w 2 > Ai \cti — oto = — 



{Ai — [Cocos(« + Ao) + CiCos(«i + Ai)]} 
" eXPl (a 2 + T CI) 

A\ + [C cos(ao + Ao) + CiCOs(ai + Ai)] > 0, 



Pr( | Im W\ | > Ai | an — a = — 



{A 2 - [C sin(a + A ) + Cisin( ai + AQ]} 

{A 2 + [Cosin(a + A ) + Cisin(ai + Ai)]} 2 



;::p| (o 2 + rcb 



A 2 ± [Cosin(a + A ) + Cisin(ai + \i)] > 0, Al + Al = A 2 , 



and 



i I iAi-Cm) 2 \ , n ( (A 2 -C M2 )'' 
Pr(|u;2|>i4) <2exp ^—"2 — + 2 ex P _ 2 , _ 2 — j» 



o' + ot J - \ o' + ot 

where 
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Cm ■ max{[C cos(ao + A ) + CiCOs(ai + Ai)]} 
a 1 -ao = ^, 

Cm = max{[Cosin(a + A ) + Cisin(«i + Ai)]} 

■n 

<Xi — OLQ = 77, 

M 
o 2 - = Z"Cl. 
For zero isi, it can be shown that optimum values of Ai and A 2 are 

(it it \ 
A 2 = A cos - + -77 , 



\4 2M 

, . /«" 77" \ 

A 2 = A sin| - + tt^ I • 



v 4 2M/ 
Even when there is ISI, we shall use these values of Ai and A 2 . 
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